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We investigate a model in which Dark Matter is stabilized by means of a Z2 parity that results from 
the same non-abelian discrete flavor symmetry which accounts for the observed pattern of neutrino 
mixing. In our A4 example the standard model is extended by three extra Higgs doublets and the 
Z2 parity emerges as a remnant of the spontaneous breaking of A4 after electroweak symmetry 
breaking. We perform an analysis of the parameter space of the model consistent with electroweak 
precision tests, collider searches and perturbativity. We determine the regions compatible with the 
observed relic dark matter density and we present prospects for detection in direct as well as indirect 
Dark Matter search experiments. 
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I. INTRODUCTION 

The existence of non-baryonic Dark Matter (DM) is well established by cosmological and astro-physical probes. 
However, despite the great experimental effort over many years, its nature still remains elusive. Elucidating the long- 
standing puzzle of the nature of dark matter constitutes one of the most important challenges of modern cosmology 
and particle physics. 

The various observations and experiments, however, constrain some of its properties p] [2]. Among the most 
important requirements a DM candidate is required to satisfy are neutrality, stability over cosmological time scales, 
and agreement with the observed relic density. While the neutrality of the particle is usually easy to accomodate in 
models, its stability in general is assumed in an ad-hoc fashion. From a particle physics point of view, the stability 
suggests the existence of a symmetry that forbids the couplings that would otherwise induce the decay. Typically, 
the most common way to stabilize the DM particle is to invoke a Z2 parity, an example of which is R parity in 
supersymmetry. 

It would certainly be more appealing to motivate such a symmetry from a top-down perspective. Different mech- 
anisms have been suggested to achieve this [3J, for instance using U(l) gauge symmetries [1HS] to get for example 
R-parity in the MSSM from a U(1)b-lU\) symmetery, global symmetries, accidental symmetries [5] or custodial 
symmetry. 

A new mechanism of stabilizing the DM has been recently proposed in Ref. [5] in which DM stability originates 
from the flavor structure of the standard model. Indeed the same discrete flavor symmetry which explains the pattern 
of neutrino mixing |10j can also stabilize the dark matter . This opens an attractive link between neutrino physics 
and DM 2 ; two sectors that show a clear need for physics beyond the Standard Model. 
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The model proposed in Ref. [5] is based on an A4 symmetry extending the Higgs sector of the SM with three scalar 
doublets. After electroweak symmetry breaking two of the scalars of the model acquire vacuum expectation values 
(vev) which spontaneously break A4 leaving a residual Z2. The lightest Zi neutral odd scalar is then automatically 
stable and will be our DM candidate. 

On the other hand, the fermionic sector is extended by four right handed neutrinos which are singlets of SU (3) x 
SU(2) x U(l). Light neutrino masses are generated via a type I see-saw mechanism [T51 - H9| . obey an inverted hierarchy 
with m^3 = and vanishing reactor neutrino angle 3 . For pioneer studies on the use of A4 for neutrino physics see 

EQ. 

We study the regions in parameter space of the model where the correct dark matter relic density is reproduced 
and the constraints from accelerators are fullfilled. We then consider the prospects for direct dark matter detection 
in underground experiments. We show that the model can potentially explain the DAMA annual modulation data 
[2"2l |2"3"] as well as the excess recently found in the COGENT experiment [23]. We show that present upper limits 
on the spin independent DM scattering cross section off nucleons can already severely constrain the parameter space 
of the model. Indirect dark matter searches through astrophysical observations are not currently probing the model 
apart from some small regions of the parameter space where the dark matter annihilation cross section is enhanced 
via a Breit-Wigner resonance. 

The paper is organized as follows: in Secjll] we present the model, in Sec. |III| the constraints from collider data 
are reviewed and in Section |TV] we study the viable regions of the parameter space. In Sec. [V] and |VT| we sketch the 
prospects for direct and indirect dark matter detection. Finally, we summarize our conclusions in Sec. |VII| 

II. THE MODEL 

We now provide a concrete realization of dark matter based upon the A4 flavor symmetry adopting the simple 
type-I seesaw framewok [L5HI9] . In order to be phenomenologically viable, the model should account not only for the 
mixing angles describing the observed pattern of neutrino oscillations [TO] but also for the two independent square 
mass splittings characterizing the "solar" and "atmospheric" sectors. These points were taken into account in the 
model proposed in [9] which we now adopt, and consider its phenomenological features in more detail. The matter 
fields are assigned to irreducible representations of the group of even permutations of four objects A4 [5J, which is 
isomorphic to the symmetry group of the tetrahedron, for more details see Appendix A. The standard model Higgs 
doublet H is assigned to a singlet representation, while the three additional Higgs doublets r\ = (771 , 772 , r?3 ) transform 
as an A4 triplet, namely 77 ~ 3. The model has in total four Higgs doublets, implying the existence of four CP even 
neutral scalars, three physical pseudoscalars, and three physical charged scalar bosons. 

In the fermion sector we have four right-handed neutrinos; three transforming as an A4 triplet Nt = (Ni, N2, N3), 
and one singlet N4. Quarks are assigned trivially as A4 singlets, as a result of which there are no predictions on their 
mixing matrix, as it is difficult to reconcile quarks with the neutrino sector 4 . 

The lepton and Higgs assignments are summarized in table]!] 
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TABLE I: Summary of the relevant quantum numbers 



3 For a similar realisation see I20| . 

4 The formulation of a full theory of flavour is beyond the scope of this work and will be left for the future. 
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The resulting leptonic Yukawa Lagrangian is: 

C = y e LJ c e H + y^L^H + y T L T l c T H + 

+y v 1 L e {N T i 1 ) 1 + y^L^N T r)) v > + y^L T {N TV ) v 
+ylL e N 4 H + M X N T N T + M 2 N 4 N A + h.c. 



(1) 



This way the held H is responsible for quark and charged lepton masses, the latter automatically diagonal. The scalar 
potential is: 

+ ^iv^v^i'ivvh" + ^'[i) t i) t ]i"[ii]i' + A5[r?V]i[w]i + ^([rfvhiirf'nh! +h.c.) 
+ Mv f vK[v f vh 2 + AgfoVkfwk + H^v]v[h^h] + \ w [rf H} 3l [H^ v ] 3l 

+ Au^VhiTi? + h.c.) + A 12 ([r 7 t ?7 t]3 i [rj H] 3i + h.c.) + A^foVkfo^k + h.c.) + X^vh^ H + h.c.) 
+ Xisdv^h^ H + h.c.) (2) 



where [..-k 2 is the product of two triplets contracted into one of the two triplet representations of A4, see eq. ( A4|,and 



is the product of two triplets contracted into a singlet representation of A4. In what follows we assume, 
for simplicity, CP conservation, so that all the couplings in the potential are real. For convenience we also assume 
A4 = A4 in order to have manifest conservation of CP in our chosen A4 basis 5 . 
The minimization of the scalar potential results in : 



(H°) = v H f 0, (r??) = v n 0, (%%) = . 



(3) 



where all vevs are real. This vev alignment breaks the group A4 to its subgroup Z 2 responsible for the stability of 
the DM as well as for the neutrino phenomenology [U] • In Appendix [B] we comment on a possible embedding of the 
model into the grand unified group SU(5). 

The stability of the DM 

Since there are no couplings with charged fermions nor quarks because of the A± symmetry, the only Yukawa 
interactions of the lightest neutral component of 772,3 are with the heavy SU(3) x SU(2) x U{1) singlet right-handed 
neutrinos. This state is charged under the Z 2 parity that survives after the spontaneous breaking of the flavor 
symmetry. One finds that, as a consequence of this Z 2 symmetry, the mass matrix for the neutral scalars is block- 



diagonal, see eq. (C8), so that the lightest neutral component of 772,3 is not mixed with the two Higgs scalars that take 



vev, H and 771. Thus quartic couplings will not induce decays for this DM candidate which is therefore stable and 
constitutes our DM candidate. 

We now show the origin of such a parity symmetry. 

As explained in Appendix |aJ the group A4, has two generators: S, and T, that satisfy the relations S 2 — T 3 — 
(ST) 3 = I. In the three dimensional basis S is given by 

/ 1 



s 



(4) 
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S is the generator of the Z 2 subgroup of A4. The alignment (77) ~ (1,0,0) breaks spontaneously A4 to Z 2 since 
(1, 0, 0) is manifestly invariant under the 5* generator: 

S(r,) = (r l ). (5) 



5 To see this, consider for instance the coupling (uiXi + \' A )(ri\ri2) 2 + h.c. arising from the terms proportional to A4 and in eq. 
Since ui + lj 2 = — 1 this coupling is real if and only if A4 = A 4 . 
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For a generic triplet irreducible representation of A4, "J = (ai, a 2 , a3) T , we have: 

/ 1 
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The Z 2 residual symmetry is defined as 
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iVi -> +Ni, 
N 2 -> —N 2 , 
N 3 -> -JV S , 



-»/2, 
-173, 



(7) 



and the rest of the matter fields are Z 2 even, because the singlet representation transforms trivially under S, see 
Appendix |Aj This is the residual symmetry which is responsible for the stability of our DM candidate. 

Neutrino phenomenology 

Here we summarize the main results obtained in Ref. [9] concerning the neutrino phenomenology. The model contains 
four heavy right-handed neutrinos so it is a special case, called (3,4), of the general type-I seesaw mechanism |19j . 
After electroweak symmetry breaking, it is characterized by Dirac and Majorana mass-matrix: 
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(8) 



where x\,X2,x$ and x 4 are respectively proportional to y", y%, y% and y\ of eq. ([!]) and are of the order of the 
electroweak scale, while Mi 2 are assumed close to the unification scale. Light neutrinos get Majorana masses by 
means of the type-I seesaw relation and the light-neutrinos mass matrix has the form: 



-m D3xi M R4> 



m n 

1 L*3 




(9) 



This texture of the light neutrino mass matrix has a null eigenvalue m.3 — corresponding to the eigenvector 
(0, —b/c, 1) T , 6 implying a vanishing reactor mixing angle #13 = and inverse hierarchy. The atmospheric an- 
gle, the solar angle and the two square mass differences can be fitted. The model implies a neutrinoless double beta 
decay effective mass parameter in the range 0.03 to 0.05 eV at 3 a, within reach of upcoming experiments. 

Notation 

After electroweak symmetry breaking and the minimization of the potential we can write: 



H' + 



H [ (v H +H' + iA' )/V2 ) ' Vl I K +H{ + iA[)/V2 J ' 



(10) 
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(H' 3 + iA' 3 )/V2 



6 Note that if we were to stick to the minimal (3,3)-type-I seesaw scheme, with just 3 SU(2) singlet states, one would find a projective 
nature of the effective tree-level light neutrino mass matrix with two zero eigenvalues, hence phenomenologically inconsistent. That is 
why we adopted the (3,4) scheme. 
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The structure of the neutral and charged scalar mass matrices follows from the exact Z 2 symmetry, which forbids 
mixings between particles with different Z 2 parities and from CP conservation. As a result the charged and neutral 
scalar mass matrices decompose into two-by-two mass matrices which, after diagonalization, give the scalar mass 
spectrum of the model. In what follows, unprimed particles denote mass eigenstates. We refer to Appendix [C] for 
details about the mass matrices and for a complete description of the mass spectrum. The Z 2 -odd sector contains 
two real CP even scalars, H 2 and H 3 , two real CP odd scalars, A 2 and A3, and four charged scalars, H 2 and H 3 . 
The Z 2 even scalars consist of two real CP even scalars H and Ho, that we generically call 'Higgses', a pseudoscalar 
Aq and two charged scalars Hq . The masses of the and Z gauge bosons impose the relation v 2 H + v 2 = v 2 where 
v stands for the standard model value of the vev. We call tan/3 the ratio of the two vevs: tan/3 = vh/v v . 

The DM candidate of the model corresponds to the lightest Z2-odd neutral spin zero particle which, for the sake 
of definiteness, we take as the CP-even state H 2 . We remind that the parameters of the model relevant for the 
DM phenomenology are the 15 couplings of the scalar potential, A^, and the ratio of the vacuum expectation values 
tan j3 = vh /v^. Indeed the minimization of the scalar potential and electroweak symmetry breaking allow to recast the 
mass parameters /ih and /i^ in terms of the couplings Ai and tan /3. Note that the couplings and the Majorana mass 
parameters in Eq. Refjl] determine the neutrino mass matrix, they are not relevant for the dark matter phenomenogy. 
Before moving to the calculation of its relic abundance in the next section we consider the phenomenological constraints 
on the parameter space of the model. 

III. PHENOMENOLOGICAL CONSTRAINTS 

In order to find the viable regions in parameter space where to perform our study of dark matter, we must impose 
the following constraints to the model : 

• Electroweak precision tests 

It is well-known that the oblique parameters S, T, U provide stringent constraints on theories beyond the Stan- 
dard Model [25] . Concerning the S and U parameters, these receive negligible contributions from the scalars 
of the model [5S1[57], hence we focus on the T parameter. We compute the effect on T induced by the scalars 
following [55] and we impose the bounds from electroweak measurements |29j : 

-0.08 <T< 0.14. 

While this bound favors a light Standard Model Higgs boson, large Higgs masses are possible in the presence 
of new physics, such as our multi-Higgs-doublet model. Indeed, a negative deviation of the electroweak T 
parameter induced by a heavy Higgs can be compensated by a positive AT produced by new scalar particles of 
the model, so one can raise the Higgs mass up to ~ 600 GeV or so (see below) 7 . 

We have explicitly verified that we can choose the mass spectrum of the model in such a way that this constraint 
is always respected. We refer to Appendix [D] for more detail. 

• Collider bounds 

Searches for supersymmetric particles at LEP place a lower bound on the chargino mass of ~ 100 GeV. To be 
conservative we apply this constraint also to the masses of the charged scalars of our model, even if slightly 
lower masses, 70 — 90 GeV, might still be consistent with LEP data|30 . 

The bounds imposed by LEP II on the masses of the neutral scalars in our model are similar to those constraining 
the Inert Doublet Model, given in Ref.[31j. This analysis applies directly to the Z 2 -odd scalar sector of our 
model of particular interest for DM phenomenology, as it constrains the mass difference between H 2 and A 2 . 
The excluded region that we adopt is taken from Fig. 7 of Rcf. 3 Lj . 



7 A similar situation holds, for instance, in the Inert Doublet Model 27 . 
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FIG. 1: Left: Tree level Feynman diagrams for the H2 annihilation into fermions. Center and right: diagrams for H2 co- 
annihilation with the pseudscalar A 2 into fermions. 

The masses of the Z 2 even neutral scalars are also constrained by LEP searches. Although a precise bound 
requires a detailed analysis of the even sector, we just impose a lower limit on these masses of 114 GeV, which 
is approximatively the LEP limit on the SM higgs mass. 

We remark that lepton flavor violating processes are suppressed by the large right-handed neutrino scale. 

• Perturbativity and vacuum stability 

The requirement of perturbativity imposes the following bounds on the Yukawa couplings of the model and 
tan/3 : 

i=l,...,15, 
tan/3 > 0.5 

this leads to an upper bound on the masses of the scalars at ~ 600 GeV 8 The constraint on tan j3 is necessary 
in order to preserve a perturbative top- Yukawa coupling. 

Finally we must impose the stability of the vacuum. The conditions ensuring that the potential is bounded from 
below are: 

Ai > , X2 + X3 + 2A 4 + A 5 > , 
Ai + 3(A 2 + A 3 + 2A 4 + A 5 ) + 3(A 9 + Mini) + 3(2A 2 - A 3 + A 8 + Min2) - 6Q > 0, 

where 

Mini = Min(X w - 2|Au|,0), 
Mini = Min(X 7 - 2|A 4 - A 5 - A 6 |, 0), 
Q = |Ai 2 | + |Ai 3 | + |A 14 | + |Ai 5 |. 

IV. RELIC DENSITY 



Note that this bound would not apply to the Inert Higgs Doublet model, which would potentially allow heavier dark matter masses. 
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FIG. 2: Tree level Feynman diagrams for the H2 annihilation into W and Z Z. 
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FIG. 3: Regions in the plane DM mass (A4dm) - lightest Higgs boson Mh allowed by collider constraints and leading to a DM 
relic abundance compatible with WMAP measurements. 



The thermal relic abundance of H2 is controlled by its annihilation cross section into SM particles. In Fig{T]and Figj2] 
we show the Feynman diagrams for the most relevant processes. In order to study the viable regions of the model we 
perform a random scan over the 16-dimcnsional parameter space (Xi and tan/3 ) and compute the dark matter relic 
abundance using the micrOMEGAs package 32, 33J. Using the mass relations in equations C12 we can trade some of 



the couplings A^ for the scalar masses. This system of 10 linear equations allows us to trade only 8 of the A^. Indeed, 
the texture of the mass matrix of the Z2 odd sector imposes the following constraints on the squared mass differences: 



ML 



3(M 
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Mjj + ) . We linearly sample on the 8 independent masses in such a way to 



fulfill the collider constraints discussed in Sec III The remaining 7 free and cos (3 are linearly sampled inside their 



allowed ranges, specified in Sec III Then, we select those choices which satisfy the perturbativity constraints for all 
the 15 Xi, the vacuum stability bounds and for which the electroweak precision tests are satisfied. Finally, we choose 
only those models which provide a relic adundance consistent with the WMAP measurements: 

0.09 < VLh 2 < 0.13. 

In Fig. [3] we show the regions with a correct relic abundance in the plane DM mass {Mum) and the lightest Higgs 
boson mass Mh ■ For dark matter masses well below the W^- threshold, dark matter annihilations into fermions are 
driven by the s-channel exchange of the Higgs scalars of the model, as shown in Fig. [TJ 

For Mh ;$ 400 GeV the annihilation cross-sections are large enough to obtain the correct relic density for all DM 
masses up to the the W ± threshold. At larger Higgs boson masses, annihilations into light fermions are suppressed so 
that the relic abundance is typically too large unless efficient co-annihilations with the pseudoscalar A 2 or with 
occur. These processes, shown in Figfljfor A2, take place only for small mass splittings between Hi and A2 or H^. 



Note that the possibility to co-annihilate with the charged scalar is ruled out, since LEP data requires Mjj > 100 
GeV. There is, however, a narrow window still allowed by LEPII limits on the A2-H2 planc[31] where arbitrarily 
small mass splitting can exist between H2 and A2 for Mom 40 GeV. For lighter dark matter particles, strong 
co-annihilations can not be reconciled with LEPII constraints as these require Ma 2 > 100 GeV [31]. We can therefore 
exclude by cosmological observations the parameter region correspondig to (simultaneously) M DM < 40 GeV and 
Mjj > 400 GeV, because the DM would be overabundant. This is clearly seen in Fig. [3] In contrast, the absence of 
points on the strip corresponding to the line Mdm ~ Mff/2 is associated to the presence of the H resonance, which 
would enhance the DM annihilation cross section giving a too small Dark Matter abundance. 

For dark matter masses larger than M\y, unless the dark matter candidate is very heavy Mdm <; 500 GeV, the 
annihilation cross section into gauge bosons is typically too large so that H 2 can only be a subdominant component of 
the dark matter budget of the Universe. However, for certain combinations of masses and parameters, the annihilations 
into gauge bosons may be suppressed by a cancellation between the Feynman diagrams (Figj2]) , leading to an acceptable 
relic density. Indeed, this happens for some points in Fig. [3] 9 . Ref. [31] shows that these processes allow for dark 
matter masses up to 160 GeV, just below the threshold for annihilations into top pairs. In our scan, the viable region 
of the parameter space extends only up to Mdm ~ 100 GeV. Solutions at higher masses cannot be excluded, though 
their scrutiny is rather involved due to the large number of parameters of the model. 

We now turn to the region Mdm ~ 500 GeV. As noted in [9], a scalar dark matter candidate annihilating into 
massive vector bosons inherits the correct relic abundance for this value of the mass if all other annihilation processes 
are absent. This scenario could be realised in our model by tuning the Yukawa couplings in order to suppress the 
dark matter annihilation into Higgs scalars and fermions. However, since the dark matter mass arises entirely from 
the electroweak symmetry breaking sector, high dark matter masses tend also to require large scalar couplings. This 
argument suggests that solutions with a heavy dark matter candidate and small couplings with the other scalars 
would be rather fine-tuned and we do not consider this possibility any further 10 . In the next two sections we study 
the prospects for direct and indirect dark matter detection. 

V. DIRECT DETECTION 



Dark Matter can be searched for in underground detectors looking for nuclear recoils induced by dark matter scattering 
against the target material. The scalar dark matter candidate we are considering couples to quarks via Higgs boson 
exchange, leading therefore to pure spin independent (SI) interactions with the nucleons. In the left panel of Fig|4]we 
show the SI scattering cross section off proton for the models with a correct dark matter abundance. We note that a 
large region of the parameter space is ruled out by the constrained imposed by current dark matter direct detection 
experiments. 

A positive signal of dark matter detection has been claimed by the DAMA collaboration [53] . DAMA has reported 
a high statistical evidence for annual modulation of the event rate over 13 year cycles [53] [53J. These results have 
prompted many attempts to interpret the data in terms of dark matter interactions with nuclei. Assuming an 
elastic WIMP interactions off nuclei and for "standard" astrophysical assumption on the local DM density and 
velocity distribution, the DAMA signal is in conflict with the null results reported by other experiments [371 l4"4l 145) . 
However astrophysical inputs are subject to large uncertainties (see e.g. |46tf48] ) and, moreover, the detector response 
is not completely known (for example the role of channelling is still under debate 49 ). We also note that there 
exist alternative particle physics scenarios where a compability between the DAMA signal and the results of other 
experiments can be obtained, e.g. [5"0"rl5"2] . 



For the Inert Doublet model, these cancellations have been noted in I31| and have been recentely studied in detail in 1341 . 
Note also that in the range ~ 60 GeV-Myy, annihilations into three body final states may play a role, as shown for the Inert Doublet 
model in I35| . 
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FIG. 4: Left plot: Spin-independent DM scattering cross section off-protons as a function of the dark matter mass. The orange 
regions delimited by the dashed (solid) line show the DAM A/LIBRA annual modulation regions including (neglecting) the 
channeling effect [36]. The green region corresponds to the COGENT data [37]. Dashed and dotted red lines correspond to the 
upper bound from CDMS (respectively from [38] and |39|1. XENON100 bounds [40] are shown as a solid red line. Right plot: 
annihilation cross section times velocity as a function of the dark matter mass. The solid and dashed red lines show respectively 
the upper bound inferred by Fermi-LAT observations of the Draco dwarf galaxy [41] and FERMFLAT measurements of the 
isotropic diffuse gamma-ray emission [42]. Projected 5 years sensitivity from measurements of the isotropic diffuse gamma-ray 
emission are shown as a dotted red line |43| . 

Recently, a hint of a possible signal from dark matter has been reported by the COGENT experiment [23]. In Figj4] 
we show the combination of SI scattering cross section and DM masses wich can fit the DAMA and COGENT results. 
An excess of events over the expected background has also been found by the CDMS experiment [38 , although with 
a low statistical significance and by CRESST [53], even if these results are still preliminary. Interestingly, all these 
possible signals point to the same region of the parameter space, favoring low dark matter masses M^m < 15 GeV 
[551 [571 HH [S3] ■ As seen in Fig|4j these possible hints of dark matter detection can be accomodated within our model. 

Current upper bounds severely challenge a possible interpretation of the aforementioned anomalies in terms of dark 
matter SI interactions. Indeed, as shown in Fig|4j the bulk of the COGENT and DAMA regions are excluded by the 
constraints inferred by XENON100 [40 and CDMS 38, 39]. Still, experimental uncertainties may significantly affect 
the upper limits obtained by direct detection experiments, expecially for low mass WIMPs (see [37j [54] [55] ) . For this 
reason a possible dark matter spin-independent interpretation of the excess, even if disfavored, cannot be excluded at 
present. 

VI. INDIRECT DETECTION 

Dark Matter indirect detection experiments look for signatures of DM annihilation into photons, neutrinos and 
(anti-) matter fluxes. The expected DM signals depend on the astrophysical details related to the DM density 
distribution in the region of observation. Particle physics enters in the determination of the DM mass, annihilation 
cross section av and the rates into various annihilation channels. In Fig. [4]we show (av) at small velocity, relevant for 
DM annihilations inside our galaxy, as a function of Mdm- For Mom ^$ 40 GeV (av) remains close to the thermal 
value at freeze-out, 3 x 10 -26 cm 3 s _1 , as expected for the s-wave DM annihilation into light fermions. At larger DM 
masses, the presence of co-annihilations allows for much smaller values of (av). The solutions at (av) ~ 10~ 25 cm 3 s _1 
correspond to DM masses just below the Higgs resonance. In this case the annihilation cross section at small velocities 
is boosted with respects to its values at the DM freeze-out. This behaviour of (av) close to a narrow Breit-Wigner 
resonance has been recently widely exploited in order to boost the annihilation signal so to explain the cosmic-rays 
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anomalies reported by the PAMELA collaboration [551159] . 

In order to sketch the prospects for indirect DM detection we show in Fig. [4] the constraints on (erf) imposed by the 
Fermi-LAT observations of the Draco dwarf spheroidal galaxy [41] and the Fermi-LAT measurements of the isotropic 
diffuse gamma-ray emission [42j . We caution that these upper bounds have been computed assuming DM annihilations 
into bb, therefore they would directly apply only for parameter choices in our model where this annihilation channel 
dominantes. Still, this happens in large part of the parameter space, and in particular at low dark matter masses, 
where the Fermi-LAT constraints are close to the predictions of the model. For a comparison of these bounds with 
those obtained for different annihiliation channels we refer the reader to the original references. Further constraints 
for different targets of observations are obtained in Ref. [60H64] 

One sees from Fig. [4] that current bounds are not yet able to significantly constrain the model. However, the Fermi- 
LAT sensitivity is expected to improve considerably with larger statistics and for different targets of observations, see 
e.g. (131 155M5] ). For example, in Fig. [4] we show the forecasted 5 years FERMI-LAT sensitivities from the isotropic 
diffuse gamma-ray emission |43j . Fermi-LAT measurements should be able to test the model for low dark matter 
masses. 

VII. CONCLUSIONS AND DISCUSSION 

We have studied a model where the stability of the dark matter particle arises from a flavor symmetry. The A4 
non-abelian discrete group accounts both for the observed pattern of neutrino mixing as well as for DM stability. We 
have analysed the constraints that follow from electroweak precision tests, collider searches and perturbativity. Relic 
dark matter density constraints exclude the region of the parameter space where simultaneously Mom ^ 40 GeV 
and Mh > 400 GeV because of the resulting over-abundance of dark matter. We have also analysed the prospects 
for direct and indirect dark matter detection and found that, although the former already excludes a large region 
in parameter space, we cannot constrain the mass of the DM candidate. In contrast, indirect DM detection is not 
yet sensitive enough to probe our predictions. However, forecasted sensitivities indicate that Fermi-LAT should start 
probing them in the near future. 

All of the above relies mainly on the properties of the scalar sector responsible for the breaking of the gauge and 
flavour symmetry. A basic idea of our approach is to link the origin of dark matter to the origin of neutrino mass 
and the understanding of the pattern of neutrino mixing, two of the most oustanding challanges in particle physics 
today. At this level one may ask what are the possible tests of this idea in the neutrino sector. Within the simplest 
scheme described in Ref. [9] one finds an inverted neutrino mass hierarchy, hence a neutrinoless double beta decay 
rate accessible to upcoming searches, while #13 = giving no CP violation in neutrino oscillations. Note however that 
the connection of dark matter to neutrino properties depends strongly on how the symmetry breaking sector couples 
to the leptons. 
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Appendix A: The A4 group 

All finite groups are completly characterized by means of a set of elements called generators of the group and a set 
of relations, so that all the elements of the group are given as product of the generators. The group A 4 consists of 
the even permutations of four objects and then contains 4!/2 = 12 elements. The generators are S and T with the 
relations S 2 — T 3 — (ST) 3 = 1, then the elements are 1, S, T, ST, TS, T 2 , ST 2 , STS, TST, T 2 S, TST 2 ,T 2 ST. 



11 





Oi = {1} 


C2 = {T} 


C 3 = {T 2 } 


C 4 = {5} 


1 


1 


1 


1 


1 


1' 


1 


UJ 


2 

W 


1 


1" 


1 




UJ 


1 


3 


3 








-1 



TABLE II: Character table of A4 where d are the different classes and uj 3 = 1. 



A4 has four irreducible representations (see Table]!!]), three singlets 1, 1' and 1" and one triplet. The one-dimensional 
unitary representations are obtained by: 



1 5 = 1 T = 1 

V S = 1 T = uj 
1" 5 = 1 T = uj 2 

where uj 3 — 1. The product rule for the singlets are: 

1 x 1 = V x 1" = 1 

V x V = 1" 

1" X 1" = 1' 

In the basis where 5 is real diagonal, 



(Al) 



(A2) 





( 1 \ 




(0 


1 




5 = 


0-10 


; t = 








!) 




i» -1/ 




\i 








(A3) 



one has the following triplet multiplication rules, 

(06)1 = 
(06)1/ 
(o6)i« 

(06)31 
(a&) 32 



= ai^i 
= ai&i 



0262 + 0363 ; 
wa 2 6 2 + w 2 a 3 6 3 ; 
u) 2 a 2 b 2 + wa 3 6 3 ; 

(a 2 63,a3&i,ai6 2 ) ; 
(a 3 6 2 ,ai63,a 2 6i) , 



(A4) 



where a = (a 1; a 2 , a 3 ) and 6 = (61, & 2 , 63). 



Appendix B: SU(5) embedding of the model 

In [5] the quark sector has not been studied and it is assumed that quarks are generically singlets of A4 in order 
to not couple to the DM. It is possible to extend such a model to the quarks by embedding it into the grand unified 
group SU(5). It is beyond our scope to give a complete gran-unified model and we only sketch a simple way to embed 
the model presented in [5] into a grand unified group. The matter assignment is the following 



Ti T 2 T 3 


F x 


F 2 


^3 


N T 


Ni 


SU(5) 10 10 10 


5 


5 


5 


1 


1 


A A 1 1' 1" 


1 


1" 


1' 


3 


1 



where we have assumed three copies of ten-multiplets and three of five-multiplets of SU(5) to decribe the three 
flavours. The scalar assignment is 

SU(5) 5 H 5 ff 5, 45 H 



A 4 113 1 
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then the Lagrangian is given by 



£down 

r 



y l { d TiF x l H + y l 2 d T 2 F 2 5 H + y l 3 4 T 3 F 3 5 H + y"' d T x Fi4£ H + y' 2 l ' d T 2 F 2 45 H + y' 3 l4 T 3 F 3 45 H ; 
V*T x T x h H + y 2 T 2 T 3 5 H + yTTxT^H + y' 2 u T 2 T 3 45 H : 

tfTiN^H + y v 2 T 2 N^ H + y^T 3 Ni5 H + y^N^ + M 1 N T N T + M 2 N 4 N 4 . 



(Bl) 
(B2) 
(B3) 



The charged lepton and down quark mass matrix are diagonal with eigenvalues 

m e = y l { d (b H ) - 3y? d (45 H ); m M = y l 2 d (b H ) - 3%' d (45ff); m T = y l 3 4 (5 H ) - 3y' 3 l ' d (45 H ); 

m d = y l i d (5 H ) +yf' d (45tf); m s = y 2 d (5 H ) +y 2 ' d (45 H )] m b = y 3 d {5 H ) + y 3 ' d {^5 H ); 



(B4) 



and the three charged lepton masses as well as the three down quark masses can be easily reproduced. The up quark 
mass matrix is 



AT 




u M ut 



M U M 



mi 
m\ 
m? 



where 



m u = y?(5 H ) ] m c = y^(5 H ) - ^ u (45 ff ); m t = y^(5 H ) + y' 2 u {45 H ). 



(B5) 



(B6) 



Given the structure of the up and down quark mass matrices, the quark mixing matrix is diagonal. While this may be 
regarded as a good first approximation, since quark mixing angles are small, clearly another ingredient is needed such 
as, possibly, radiative corrections or extra vectorlike quark states. A full fit of the quark sector observables within a 
unified extension incorporating the flavour symmetry is beyond the scope of our paper. 

Appendix C: Mass spectrum 



To simplify the notation we define the following combinations of couplings: 

L = Ag + Aio + 2An, 

Q = A12 + A13 + A14 + A15, 

P = A2 + A3 + 2A4 + A5, 

Ri = (-3A 3 -6A 4 + 2A 6 + A 7 + A 8 ), 

R2 = (-3A 3 -2A4-2A 6 + A 7 + A 8 -4A 5 ), 

Rs = (-3A 3 - 4A 4 - 2A 5 + A 8 ). 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 
(C6) 



The neutral scalars mass-matrix in the basis H' — H[ — A' — A[ — H 2 — H 3 — A' 2 — A' 3 is block diagonal because of 
the Z 2 symmetry and CP-conservation. It is given by: 



M neu trals — 



V 



HH 











M GAo 











M H2 h 3 
















Ma 2 a 3 ) 



(C7) 



and the charged scalars mass matrix in the basis 7J — — H 2 — H 3 is: 
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/ M r+H + \ 

M eharged =[ G+H ° . (C8) 

The matrices M G A a an d M G+H + have a vanishing eigenvalue corresponding respectively to the neutral and charged 
Goldstone bosons eaten up by the Z and W ± gauge bosons. The diagonalization of the mass matrices goes as follows: 

Dhh = uIj^Mhh^hh^ D GAq = U GAq M GAo Uga ] 

D G+H + = U G+H +M G+H +U G+H +; M X2X3 = tf 23 M X2X3 U 23 , (C9) 

with X = H, A, H+ and 

UGA ° = Ug+h * ~ [ sin/3 cos/3 ) ' U "-{ 1/V2 1/V2 ) ' (C10) 

The mixing angle between 772 and 773 is 7r/4. The ratio of the vevs was parametrized by the angle j3 that controls the 

mixing within the z^-evcn sector for the charged scalars and pseudoscalars: 

2 2 , 2 
v =v v + v H ; 

tan (3 = vh /v v . (Cll) 
We do not report the lenghty expression for Uhh - The mass spectrum reads then as: 

M\ a = -2w 2 An, (C12) 

Af| = XiVh 2 + Pv 2 - ^f{v 2 v H 2 {L 2 - 4PA0 + (Pv, 2 + A^ 2 ) 2 ), (C13) 

M 2 Ha = X lVH 2 + Pv v 2 + y/{v n 2 v H 2 (L 2 - APX^ + (Pv v 2 + X lVH 2 ) 2 ), (C14) 

M 2 + = -(Aio + 2A 11 )u 2 /2 ! (C15) 

Mf 2 = v ri (-3Qv H + R lVv )/2, (C16) 

M 2 2 = (-4A n V + Qv v v H + R 2 v n 2 )/2, (C17) 

M 2 = (v H 2 (X w + 2X^) - Qv v v H + R 3 v v 2 )/2, (C18) 

M| 3 = M| 2 +3guH Uj7 , (C19) 

M 2 3 = Af 2 2 + Q^, (C20) 

M i+ = M l+ + Q v hVv ( C21 ) 

Appendix D: Oblique parameters 

Following the notation of [28], the T oblique parameter for the standard model extended by n higgs doublets with 
hypercharge 1/2 is : 



2 ( n 2n 

117 La=2 6=2 

2n-l 2n 

"E E [lm(V^)J 2 F( M 2 , M 2 ) (D2) 

Ti—l n 

-2E E K^U^K^') (D3) 

a=2 a'— a+1 
2n 

+3 £ [Im J 2 [F (m|, /^ 2 ) - F (m 2 ,, /J 2 )] (D4) 

6=2 

-3[f(m|,m 2 ft )^(m^mJ)]}, (D5) 
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where m a , m a denote the masses of the charged scalars and /i^, are the masses of the neutral ones, a is the 
fine-structure constant and the function F is defined as (x,y > 0): 

!x + y xy x 
— x In - <*= x * y, 
2 x-y y (D6) 
<= x = y. 

We evaluate the U and V matrices for our model as: 



U = U <**t ° • V= iUc+H t UHH ° ° ° (D7) 

V o u 23 ) ' 1 o o u 23 iu 23 J 
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